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1. Introduction 

As the precision of gravitational experiments increases it is increasingly necessary to 
compute with high accuracy the gravitational fields of complicated objects. This paper 
provides analytic results that have proved useful for this task. 

Multipole expansions provide a powerful formalism useful in many physical 
applications. One such application occurs in equivalence principle tests P P where 
gravity gradients set a fundamental limit to the precision of the results. Spurious 
effects from gravity gradients can be minimized in a systematic way by nulling the 
low-order multipole moments of the differential accelerometer (typically a torsion 
pendulum) and the attractor. In the case of torsion balances, where differential 
acceleration of the test bodies produces a torque on the pendulum, expansions in 
spherical multipoles are particularly useful. Evaluating the multipole moments of 
complicated objects presents a challenge. The moments can be found by decomposing 
the objects into elementary 3-dimensional shapes whose moments can be calculated 
exactly. Because the moments have definite rotational and translational properties^, 
a catalog of moments of the elementary solids evaluated in any suitable coordinate 
system allows one to “assemble” realistic objects and compute their multipole moments 
with high precision. The formalism of Ref. 0, together with the moments cataloged 
here, provides a convenient way to understand the gravitational response of torsion 
balances, including the effects of imperfections in the instruments themselves. This 
paper gives analytic expressions for the inner multipole moments of 11 elementary 
solid objects that have proved useful in modeling the gravitational response of actual 
devices. 

Another application of analytic results is found in the gravitational compensation 
of drag-free spacecraft. For example, the LISA (Laser Interferometer Space Antenna) 
gravitational wave detector will use three such spacecraft, each containing two inertial 
bodies that form the end mirrors of interferometer arms. Stringent requirements on the 
test-body DC acceleration and force noise require accurate modeling of gravitational 
forces, torques, and gradients produced by the surrounding spacecraft enclosure |lj. 
Since much of the disturbing mass distribution is in close proximity to the inertial 
bodies, the multipole formalism is not valid for this task. The analytic expression for 
the force between two rectangular solids, derived here, provides an accurate tool for 
assessing these gravitational disturbances. 

2. Inner moments of the elementary solids 

We present analytic expressions for the gravitational (or Coulomb) inner moments, 
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of the following elementary shapes: sphere, rectangular solid, cylindrical annulus, 
cone, triangular slab, trapezoidal slab, cylindrical section, right tetrahedron, hole in 
a cylinder, hole in a plate, and pyramid. All moments are given with respect to the 
origin, and are proportional to the mass density, p. We list all non-vanishing moments 
with l < 5; the unlisted l < 5 moments vanish identically. 
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2.1. Sphere 

The sphere of radius r is centered on the origin. 

2 

<7oo = -Px/tt?’ 3 

2.2. Rectangular box 

The box is centered on the origin and is specified by the positive quantities x, y, and 
2 !—its dimensions along x, y, and z. 
pxyz 

900 = 

<720 = -900 ( X 2 + y 2 - 2 z 2 ) 

<722 =900^^1 (x 2 - y 2 ) 

940 = 9 oo ^ [9 .t 4 + 10 x 2 y 2 + 9 y 4 - 40 (ar 2 + y 2 ) z 2 + 24z 4 ] 

942 = ( x2 - y ' 2 ) [ 3 ( x2 +^ 2 ) - 10z2 ] 

944 = 900 1^8 ^ ~ 10x2y2 + 3y4 ) 

2.3. Cylindrical annulus 

The cylindrical annulus has its axis of symmetry along z and is symmetric under 
reflections through the x-y plane. It is specified by inner and outer radii, r i and r 2 , 
and height h > 0. It need not be a complete annulus, but spans azimuthal angles 
between (fix and <fi 2 , which are measured with respect to x. We clehne A$ = (fi 2 — (fii, 
4> = {(fii+ <fi 2 )/ 2 , A R = r 2 - n, and r = (n + r 2 )/2. 

9oo = —^ftrA$A fl 

zy 7T 

9n = ~J^ e ~ l<PhAR + 12r2 ) sin 

920 = ~^\[~ hr A-J-A/j (3A| - 2ft, 2 + 12r 2 ) 

922 = ^77 \f^-e~ 2l(t, hrA R (A r + 4r 2 ) sin A$ 
lo V Z7r 

9 31 = \f^ 6 ~ l(l>hAR 9A ^ “ 20A ^ ^ 2 “ l8r2 ) _ 240? ’ 2 _ 3r2 ) Sh4 

933 = “i k\fi e ~ 3i ' hAR ^ R+40A « r2 +80r4) sin ^ 

9 40 = /— h r A<x,A R 15A|j — 40Afj (ft 2 — 5?’ 2 ) + 8 (ft 4 — 20 h 2 r 2 + 30r 4 ) 

lzo(Jy7r L 

942 = “ 128 \/ S e ” 2i ^ rAii [ 3A ^ " ^ ^ + Ar ('~ 6h2 + 4 ° r2 ^ 


sin A$ 
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n / Q K 

944 = 3T— \ —e~ 4q, ^hrA R (2>A 4 R + 40A 2 R r 2 + 48r 4 ) sin2A$ 
O-LZ V ^7T 


951 — 


953 = 


955 = 


J^e-^hAn 


3584 V 30tt 

+56A|j (d 4 - 60 h 2 r 2 + 150r 4 ) 

P ^ 


15A| - 84A|j ( h 2 - 15r 2 ) + 672r 2 (d 4 - 10 h 2 r 2 + 10r 4 ) 
A$ 


sin ■ 


9216 V 35?r 


21 e - 3 ^/iA H 


15A| - 56A|j/i 2 + 140A| (9A| - 16d 2 ) r 2 
3A$ 


+560 (l5A^ - 8 h 2 ) r 4 + 6720 r 6 


sm ■ 


3 p 

5120 V 7tt' 


- 5i ^hA R (^A 6 R + 84A|jr 2 + 560A|r 4 + 448r 6 ) sin 


5A$ 


2-4- Cone 


The (truncated) cone is symmetrical about z. and is specified by upper and lower 
radii, r± and r 2 , and height h > 0. The cone extends a distance d/2 above and below 
the x-y plane. Complete cones have vanishing values of r± or r 2 ■ 

9oo = ^'/* h ( r i + r i r 2 + r t) 

9io = | ^h 2 (?{ + 2rir 2 + 3r 2 ) 

2 h 2 (r 2 + 3rir 2 + 6r 2 ) - 3(r 4 + r\r 2 + r 2 r| + rir 2 + r 2 ) 


920 = 


24 

P 


93° = ;^\/77rd 2 2 ^ 2 ( r i + 4rir 2 + 10r 2 ) - 3(r 4 + 2 r\r 2 + 3r 2 r 2 + 4rir 2 + 5r 2 ) 

940 = -p-\/7rd 8 h 4 (r 2 + 5rir 2 + 15r 2 ) - 12d 2 (r 4 + 3rfr 2 + 6r 2 r| + lOrirf + 15r 2 ) 

ODU 

+15 (rf + r\r 2 + r 4 r| + r 3 r 2 + r 2 r 2 + nr 2 + rf) 

950 = — ^—Vllnh 2 8d 4 (r 2 + 6r!r 2 + 21r|) - 12d 2 (r 4 + 4r 3 r 2 + 10r 2 r| + 20r 1 r 2 + 35r 2 ) 

zboo 

+15(rf + 2 r\r 2 + 3r 4 r 2 + 4r 3 r 2 + 5r 2 r 2 + 6rir 2 + 7r 2 ) 


2.5. Triangular slab 

This shape has reflection symmetry about the x-y plane and thickness h > 0. The 
triangular faces have vertices at ( x,y ) = (0,0), (d,y i) and ( d,y 2 ). The restriction 
that one side be parallel to ij is easily overcome using the rotational properties of the 
moments [3). 

9oo = ~T~/=h\d (yi - y 2 ) \ 

4+7T 

9n = 9oo [-2 d + i{yi+ y 2 )] 

920 = — 9 oo~j^- ( 3 d 2 — h~ + y\ + y\yi + y 2 ) 

l 

922 — 900 ^ 


3d 2 -y\- yiV 2 -y\- 3 id (yi + y 2 ) 
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931 — 9oo 


1 ^ 0 \j\ {36d 3 - 18 id 2 (yi + y 2 ) + i ( y\ + y 2 ) [10ft 2 -9 (yf + yf)] 
-4 d [5ft 2 -3 (yf + yry 2 + yf)] j 


1 7 


933 = 900 g y - [(1 + i) d + 2/1 - iy 2 \ x [(1 + i) d - iyi + y 2 ] [2 id + {y i + y 2 )] 

9-4° = 9oo g'g {30d 4 + 3 ft 4 - 10ft 2 (yf + y 1 y 2 + yf) + 6 (yf + yfy 2 + yfyf + V\y\ + y \) 
+10d 2 [—3ft. 2 + 2 (yf + yiy 2 + yf)] } 

942 — — 9o° {^O^ 4 — 15d 2 ft 2 — 20id 3 (3/1 + y 2 ) + 5ft 2 (yf + yiy 2 + yf) 

~4 (yf + yfy 2 + yfyf + yiyf + y 2 ) + 5 id (yi + y 2 ) [3ft 2 — 2 (yf + yf)] | 
1 / 7 r 

944 = 9oo g y |5d 4 + vi + yf y 2 + yf yl + yiyf + vi ~ !0*d 3 (yi + 2/2) 


+5*d (yi + y 2 ) (yf + vi) - 10 d 2 (yf + yry 2 + yf) 


951 — —9oo 


1 

112 V 30 


|60<i 5 - 30*ft 4 (yi + y 2 ) + 6id 2 (yi + y 2 ) [7ft 2 - 5(yf + yf)] 

-4d 3 [21ft 2 - 10 (yf + yiy 2 + yf)] + 2d[7ft 4 - 14ft 2 (yf + yiy 2 + yf) 
+6 (yf + y?y 2 + yfyl + yiyf + yf) ] 

-*[7ft 4 (yi + y 2 ) - 21ft 2 (yf + yfy 2 + yiyf + yf) 

+10 (yf + yfy 2 + yfyf + yfyf + yiyf + yf) 


953 = 900 48 Y 35 { 3 °d 5 _ 45 * rf4 + 2/ 2 ) + 3 *d 2 (J / 1 + 2 / 2 ) [l4ft 2 - 5 (yf + ; 

-4d 3 [7ft 2 + 5 (yf + yiy 2 + yf)] 

-* (yi + 2 / 2 ) [7ft 2 (yf + yl) - 5 (yf + yfyf + yf)] 

+d 28ft 2 (yf + yiy 2 + yf) - 18 (yf + yfy 2 + yfyf + yiyf + yf) j 

1 fir 


955 — 900 


16 V 7 


[—2d + i (yi + y 2 )] [d 2 - yf + yiy 2 - yl - id (yi + y 2 )] 
x [3d 2 - yf - yiy 2 - yl - 3fd (yi + y 2 )] 


ft. Trapezoidal slab 

The trapezoidal slab has reflection symmetry about the x-y plane and a thickness 
t > 0. Two of its sides are parallel to y with lengths Wi and w 2 . It is symmetric about 
x. i.e., the two parallel sides are centered at y = 0 with separation ft > 0. We define 
w = (w 1 + w 2 )/ 2 and A w = (w\ — w 2 )/w. 

9 


9oo = 


2 y/n 


htw 


1 


9n = q 00 ^h(A w - 6) 
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<720 = 9 oo ^ [4A w h 2 - A lw 2 - 4 (4ft 2 - 2 1 2 + in 2 )] 

922 = ~ 9 oo^y^ [4A ffi /j 2 + A^w 2 + 4 (—4ft 2 + w 2 )] 

ffii = -goo^ [3Aj > 2 - 30A> 2 - 120 ( 6 ft 2 - 4t 2 + in 2 ) 

+4A„, (54ft 2 - 20t 2 + 15 m 2 ) 

933 = 9oo ~~ A ™ w2 + 10A ^ w2 + 4 A W ( 6 ft 2 - 5in 2 ) +40 (-2ft 2 + in 2 ) 

g 4 o = 9 ooy^{ - 24A lw 2 h 2 + 9A* in 4 - 32A m ft 2 (24ft 2 - 20f 2 + 15in 2 ) 

+ 8 A 2 ,in 2 (24ft 2 - 20t 2 + 15in 2 ) + 16[l44ft 4 + 24f 4 - 40f 2 in 2 + 9in 4 
—40ft 2 (4i 2 - w 2 ) ] } 

942 = 9 °° 1024\/l0 32A “’ ~~ + 3A “ w4 

+40 A 2 ,in 2 (- 1 2 + in 2 ) + 16 (-4ft 2 + in 2 ) ( 12 ft 2 - 107 2 + 3 w 2 ) 

944 = 900777777 : 1/777 24A[[ in 2 ft 2 + 3A 4 in 4 + A w (-256ft 4 + 480h 2 w 2 ) 

+A [> 2 (-192ft 2 + 40in 2 ) + 16 (48ft 4 - 40ft 2 in 2 + 3w 4 ) 

95 i = q 00 -^^\j^h{3A 5 w w 4 - 42A 4 in 4 + 8 A^ ■in 2 (26ft 2 - 14t 2 

+21w 2 ) - 224 [80ft 4 + 24t 4 - 20£ 2 w 2 + 3in 4 

-20ft 2 ( 6 t 2 - w 2 ) ] - l\2A 2 w w 2 [14ft 2 + 5 (- 21 2 + w 2 )] 

+16A„, [400ft 4 - 252ft 2 (2t 2 - w 2 ) + 7 ( 8 t 4 - 20t 2 in 2 + 5in 4 ) ] | 

953 = “ 9 oo 49^52 “ 9A ^> 4 + 126A 4 w 4 + H2A 2 w 2 (14ft 2 

—20t 2 + 15in 2 ) - SA 3 w w 2 (26ft 2 - 28f 2 + 63in 2 ) 

-224 [80ft 4 + 40f 2 w 2 - 9 m 4 - 20ft 2 (4t 2 + to 2 )] 

+16A W [400ft 4 - 84ft 2 (4t 2 + 3in 2 ) + 35 ( 8 f 2 in 2 - 3w 4 )] } 

955 = 900 1 63 84 \/y 1h 3A » w4 “ 42 A t ^ 4 + A^w 2 (-208ft 2 + 168in 2 ) 

+112A 2 in 2 (14ft 2 - 5tn 2 ) + 16A W (80ft 4 - 252ft 2 tn 2 + 35in 4 ) 

-224 (16ft 4 - 20ft 2 in 2 + 3w 4 ) 

2.7. Cylindrical section 

The cylindrical section has a cross section consisting of the space between a chord 
of a circle (parallel to the y axis) and the arc the chord subtends. This object has 
reflection symmetry about the x-y plane and is specified by its thickness, t > 0 , the 
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circle radius, r, and the distance along the x axis, d > 0 , from the center of the circle 
to the chord. We define s = \J 1 — d 2 /r 2 . 

P r _i fd\ 

q°° ~ ~%Jr^r ~ sd + rC0S 

5n = ~i= Str (d 2 - r 2 ) 

V D7T 


P /5 


24 V tt 


tr Sd (2d 2 -t 2 +r 2 )+r (t 2 - 3r 2 ) cos 1 (^j 


q22= 2\h Sdtr 

q ^ = i^ Str ^- d2 ) (M 2 -5h 2 + 9r 2 ) 


q 33 = try —S ( 6 d 4 - 7d 2 r 2 + r 4 ) 

q i0 = i 6 ( |^_ tr| - Sd(l6d 4 - 20d 2 f 2 + 3t 4 + 2 (4d 2 - 5f 2 ) r 2 + 6 r 4 ) 


+3r (t 4 - 10 t 2 r 2 + 10r 4 ) cos -1 ^ | 


^42 = —- Sdtr (d 2 — r 2 ) (Ad 2 — 5 1 2 + 6 r 2 ) 

8v 107T 

<744 = ——dtr\l-J—S ( 8 d 4 — lld 2 r 2 + 3r 4 ) 

8 V 10-7T v ' 


q 51 = -j~ tr ( d2 ~ r 2 ) \E-s\wd 4 - 28 d 2 t 2 + 7t 4 + 6 (4d 2 - 7t 2 ) r 2 + 30r 4 

1IZ V OU7T L 

953 = — tr (d 2 - r 2 ) \[As\- 48 d 4 - Ut 2 r 2 + 15r 4 + 12d 2 ( 7t 2 - 6r 2 ) 
144 V 3o7r L J 


<755 = tr (d 2 - r 2 ) J—S (80 d 4 - 48 d 2 r 2 + 3r 4 ) 


£.(§. Right tetrahedron 

This shape consists of a tetrahedron having three mutually perpendicular triangular 
faces that meet at the origin. The fourth triangular face is defined by points at 
coordinates x, y and z along the x, y, and z axes, respectively. 


9io — 9oo -^-z 

du = ( x _ i V ) 

920 = ~ qoo f^ ( x 2 + y 2 - 2 z 2 ) 
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<721 =-q 00 j^z (x- iy) 

3 

922 = 900 ivM ~ ixy ~ v > 

<730 = “ 900-^2 ( X 2 + y 2 - 2z 2 ) 

1 [7 


931 [3a: 3 - ia: 2 y + a;y 2 - 3iy 3 - 4 (a; - iy) z 2 ] 

i r~7~ 

932 = 900 g y ^2 (a; 2 - ixy - y 2 ) 


933 = —900 


(x 2 - y 2 ) (x - iy) 


j_ n 

16 


9 40 = 9007^ [3x 4 + a: 2 y 2 + 3y 4 - 4 (a; 2 + y 2 ) z 2 + 8 z 4 ] 

3 

941 = 900 irvi 2 [ 3a;3 ~ * x22/ ~ 3 * y3 + Aiyz 2 +^ ~ 4z2 )] 

942 = —go° 56v /yq l 2x4 ~ ix3 y - ix y 3 - 2 v 4 ~ 2 l ^ 2 - ix y ~ y 2 ) ^ 2 ] 

3 

943 = -9 o ° i6v /^ 2 (a; 2 - y 2 ) (x - iy) 

944 = 9o° “7^ (a; 4 - ia; 3 y - x 2 y 2 + ixy 3 + y 4 ) 

950 = 90 °4^6 v/ ^‘ 2 [ 3a;4 + + 32/4 ~ 4 ^ 2 + y2 ^ 22 + ^ 


/II 


95i — — 9oo 


5a; 5 — ix 4 y + x 3 y 2 — ix 2 y 3 + xy 4 — 5 iy 5 — 


448 V 30 L 

2 (a: + iy) (3a; 2 — 4ia;y — 3y 2 ) z 2 + 8 (x — iy) 2 4 
952 = 9oo ^ y^ys [^ 2a:4 + ix3 V + + 2y 4 + 2 (a; 2 - ixy - y 2 ) z 2 ] 


1 11 


5a; 5 — 3ia; 4 y — x 3 y 2 — ix 2 y 3 — 3 xy 4 + 5iy 5 — 4 (a; 2 — y 2 ) {x — iy) z z 


953 = 900 1 28 V 35 L 

3 /II / 4 .3 22 - 3 41 2 

954 = 900gj W ^2 (a: - ix y - x y + ixy +y ) z 

3 /II / . \ ( 9 2\2 2 2 

955 = -^0° 128 V y ( x _ iy > [I 21 “SM + *2/ 


2.P. Cylinder hole 

This shape consists of the volume that would be removed by drilling a hole of radius r 
into a cylinder of radius R. The symmetry axis of the hole is along z, and the cylinder 
has its symmetry axis along y. The moments require the hypergeometric function 
2 Fi(a,b; c;x). 

/II r 2 \ 

9oo = p\f^r 2 R 2 F 1 ( --,-;2; ^2 ) 
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920 


= -^r 2 R\(r 2 -2R 2 ) 2 Fi(-± \,2-^ 


z F / _ X 3 Q 7 

2jF| 1 2 ’ 2 ’ 3 ’ _R 2 


15-7T 


922 - p\l ——r R 2F1 ( 2; — 


1 1 


940 = 40^ ri? ^ “ 4 ° r3R2 + 24rR ^ 2Fl ( -o’ 7 ^ m 


1 1 


2’2’ ’ i? 2 


1 3 


+ (l3r 5 - 32r 3 i? 2 ) 2 ^ ( 3; -^ ) + 16r 5 2 F, ( 4; — 


P / 7T 4n 
942 = 77 \ / T7T r ^ 


8 V 10 


(6r 2 — 20i? 2 ) 2 Fi(-i 1;2;-^ 


1 3 


1 5 


2’2’ ’ i? 2 


+2 (r 2 + 10i? 2 ) 2 F 1 ( -± ^; 3; ^ ) - 13r 2 2 F, ( ^;4; ^ 


1 5 


2’ 2’ ’ i? 2 


944 = 


3/3 /357T 2 


2 Fi - 


1 I 9 _ 

2 ’ 2 ’ ’i? 2 


2.70. Plate hole 

This shape consists of the volume that would be removed by drilling a hole of radius 
r through a parallel-sided plate of thickness t. The plate is centered on the x-y plane. 
The hole axis, which passes through the origin, lies in the x-z plane at an angle 6 
measured from z, where —7r/2 < 6 < 7r/2. We define s = tan61 

ptr 2 sjTt(\ + s 2 ) 


9oo = 


75 


920 = 9007 ^ [3r 2 (s 2 + 2) + 7 2 (s 2 - 2)] 


921 — — 9oo 


1 / 5 


4 V 6 

922 = - 9oo^y^s 2 (3r 2 + 6 2 ) 


t s 


940 = - qoo — [I0r 4 (8 + 8s 2 + 3s 4 ) + 10tV(3s 4 - 8) + 6 4 (8 - 24s 2 + 3s 4 )] 

94 1 = - 9oo^^6 2 s [5r 2 (4 + 3s 2 ) + t 2 (3s 2 - 4)] 


942 — 9oo 

943 = 9oo 


1 

64 V 10 
3 

64 


s 2 [lOr 4 (s 2 + 2) + 10r 2 s 2 6 2 + 1 4 (s 2 - 6)] 
f-sV (5,-» + C) 


944 = ~ j^s 4 (lOr 4 + I 0 r 2 t 2 +t 4 ) , 
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This shape consists of a symmetric pyramid whose base with side lengths x and y lies 
in the x-y plane centered about the origin, and whose apex is at ( x,y,z ) = (0,0, h). 
phxy 


9oo = 


9io — 9oo - 


6 y/TT 

VS h 


4 h? — x 2 — y 2 
8x/5 

3 x 2 — y 2 
10 8 

930 = doo^h(4:h 2 - x 2 - y 2 ) 


920 — 9oo - 


922 — 9oo 1 


932 — 900 1 


7 h(x 2 — y 2 ) 
30 16 


940 = 900 4 4 80 [ 128/l4 + 9a;4 + 10x2 y 2 + 9y 4 - S2h 2 {x 2 + y 2 )\ 

942 = q °° AAo ~ y 2 } ^ 3 ( a;2 + y 2 )\ 


944 = 9oo 


448-\/T0 
3 

128V70 

v'TI, 


[3cc 4 - 10 x 2 y 2 + 3 y 2 ] 


950 = ®o^ggft [128h 4 + (x 4 + 10:r 2 y 2 + 9y 2 - 32 h 2 (x 2 + y 2 )] 


952 — 900 


1 / 33 

512 V 70 


h{x 2 - y 2 ) [8h 2 - 3(ir 2 + y 2 )] 


954 = 900 


1024 V 70 


(3s 4 — 10s 2 y 2 + 3y 4 ) , 


3. Gravitational force between two rectangular solids at arbitrary 
separations 

In situations where the test masses are in close proximity, the multipole formulation is 
not valid. In such cases one can use a rectangular mesh to model complex objects. We 
present here an analytic expression for the force between two parallel-sided rectangular 
solids at arbitrary separations. Mass 2 has density p 2 and is centered at the origin with 
side half-lengths (a 2 ,5 2 ,c 2 ). Mass 1 has density p 2 and is centered at (Ax, Ay, As) 
with side half-lengths (ai,6i,ci). 

f z = -G PiP2 [il - Io 2 - if + if - yo (if -lo 1 ) + (yi - yo)(i 2 1 -if) + ^(if-if)], (i) 
where I 4 , if, and the y, are derived in the Appendix. Similar expressions yield F x 
and F y . Although the complete expressions are quite long they can be evaluated very 
quickly on modern computers. The precision attainable with our expressions may be 
very useful when analyzing the gravitational balancing of drag-free spacecraft, such as 
those envisioned for LISAgj. Our expression may also find use in accurately modeling 
the Newtonian forces between rectangular test masses in future short-range tests of 
gravity. 
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The z-component of the gravitational force on mass 2 is given by the 6-dimensional 
integral 


F z — —Gpip2 j 


pAx-\-a\ pAy-\-bi rAz-\-ci raz pb 2 rC 2 


d U 


/ Ax—ai J Ay—bi J Az—c\ J — a^, • 


Lc 2 d (“l - 22) 


dz 2 Ay 2 dx 2 dzi dyi dxi , 


_ 1 I c\ 

where U = [(xi — X2 ) 2 + {yi — V2 ) 2 + (21 — Z2) 2 ] . We transform the coordinates 

into sums and differences. For example, 


2 = Z\ — Z2 


Z s = 2i + Z 2 


d 2 d 2 s = 2dzi dz 2 ■ 

Integration over z s is trivial because the function is independent of this variable. The 
z-integrations give 


d2i / dZ2 


dU 


1 


dU 


d( 2 i - z 2 ) 2 


= - dz dz* = 


dz 


Udz - 


Udz 


= F 0 (x, y) + F 3 {x,y)~ Fi (x, y) - F 2 (x,y) , 

where 

Fi(x, y) = In ^ Zi + sj x 2 + y 2 + zf 

The limits are given by 

z 0 = Az — ci — c 2 

Zi = Min (Az + Ci — c 2 , Az — Ci + c 2 ) 
z 2 = Max (Az + Ci — c 2 , Az — ci + c 2 ) 

23 = Ax + ci + c 2 . 



Analogous expressions define the limits for the other coordinates. The x integrand 
involves terms such as Fi(x,y) and xFi(x,y), which lead to 


Gy (v)= / dx Fi (x, y) = xj [In (zi + ry) - 1 ]+y 


tan 


— tan 


V r ij 


+Zi In (xj + ry) 


and 


Gy (V) = J dx x Fi (x, y) = - [2 Zi n j - rjj + 2 p 2 In (z { + ry)] , 

where ry = yjx 2 + y 2 + zf and pj = yjx 2 + y 2 . After performing the trivial integrals 
over x s and y Sl five of the six integrations are complete. The final step is the integration 
over y 1 which contains the terms 


J G\(y')dy'=x- } y 


In (zi + ry) - - 


+ yzi [In (xj Try)- 1] 





Xj Zi \ 

ynj 


z- 

~i 

2 


2 tan -1 1 — 

2i 
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— tan 


-i Zj y 


Z\ r“i- 


tan 


-l y 


x\ 


-i y z i 


— tan 

.:/ V x i r ij 


+£j Zj In (y + ry) , 


J y'Glj(y')dy' = ^ {-3r? (ay + Zi) - 4ay(y 2 - ziry) + 6y 2 Zi \n( Xi + ry) 

+4y 3 [tan^f^ - tan" 1 )^) +6x j2 / 2 ln^+ry) 

36 (xf + zf) 


—2 zf In 


\ 3 | ^6(ay + z f) 

xf zf ( y 2 + zf Kay + ry) J 2 ln 1 x f zf(ry + zf) 


/ Gl{y')&y'= -^6 


tan _1 ( — tan -1 ^ 




+ Zi (3®? + zf) ln(y + ry) 


+y (3a; 2 + y 2 ) ln(zi + ry) > - y (5ry + 16a; 2 + 4zf - 12ryZi) 


and 

rv 


V ' Gfj (y') dy ' = ^ 8zf ry - (3p? + zf) 2 - 6p 2 zf + 12pf ln(z ; + ry) 

Finally, we define 

ff k (y) = G 2 ,(y) - G 2 0 (y) - a; 0 (G^y) - G 3 0 (y)) + (*i - x 0 ) (G 3 2 (y) - G^y)) 


+^(G 3 3 (y) - G^y)) - G 2 3 (y) + G 2 2 (y) 


and 


ry n 

T n= J dy [H 0 (y) + H 3 (y) - H 3 {y) - H 2 {y)\ 


ry n 


/ 2 = 


ydy [i? 0 (y) + ff 3 (y) - i?i(y) - ff 2 (y)] 


which lead to Equation ^ 
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